AD-A096  6t>5  WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH  CENTER  F/G  12/1 

ON  THE  CHOICE  OF  THE  EXTERIOR  KNOTS  IN  THE  B-SPLINE  BASIS  FOR  A— ETC(U) 
DEC  80  J  KOZAK  DAAG29-80-C-0041 

UNCLASSIFIED  MRC-TSR-2148  NL 


END 

4 -S  I 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 

December  1930 

Received  May  29,  1980 


MAR  2  3  198] 


Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 

U,  S.  Army  Research  Office 
p.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


Council  for  International 
Exchange  of  Scholars 
Suite  300,  Eleven  DuPont  Circle 
Washington,  D.  C.  20036 

81  3  19  085 


The  B-spline  representation  of  a  spline  on  some  interval  [a,b]  requires 
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linear  interpolation  conditions,  and  so  is  preferable  when  the  B-spline 
representation  of  a  spline  interpolant  is  to  be  constructed.  Such  a  choice 
usually  also  simplifies  the  algorithms.  In  particular,  one  is  able  to  compute 
stably  the  B-spline  coefficients  of  a  complete  spline  interpolant  by  Gauss 
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SIGNIFICANCE  AND  EXPLANATION 


The  B-spline  representation  of  a  spline  on  some  interval  [a,b]  requires 
the  introduction  of  additional  knots  outside  ]a,b[  that  have  nothing  in 
common  with  the  spline  space  itself.  Nevertheless,  their  choice  influences 
the  numerical  accuracy  when  the  B-spline  representation 

5  “i  'W^V 

of  a  spline  approximant  is  to  be  constructed.  In  this  note  we  prove  that 
locating  all  the  additional  knots  at  the  endpoints  a  and  b  of  the  interval 


minimizes  the  norm 


li  "i.k.t’'1'— Cl/ 


and  is  therefore  preferable.  Here  X  :=  (X  )  is  an  arbitrary  sequence  of 

linear  functionals  with  support  in  I  except  for  the  assumption  that  the 

interpolation  problem  is  correct,  i.e.  (X.  N.  )  nonsingular.  The  above 

i  3  i  X ,  t 

choice  usually  also  simplifies  the  algorithms.  In  particular,  one  is  able  to 
compute  in  a  stable  way  the  B-spline  coefficients  of  a  complete  spline 
interpolant  by  Gauss  elimination  without  pivoting  though  the  corresponding 
matrix  fails  to  be  totally  positive  or  diagonally  dominant. 


The  proof  relies  on  the  observation  that 

‘"l.*,,1  *  1Ni.k,t,0t  T  °n  Ia'M  ’ 
with  being  two  knot  sequences  that  coincide  in  [a,b].  The  matrix 

^  for  the  particular  choice  of  additional  knots  t_  described  above  turns 

to  be  totally  positive  and  its  columns  sum  up  to  1  ,  i.e.  Q  ^  is  composed 


of  (discrete)  B-spllnes. 


-X 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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ON  THE  CHOICE  OF  THE  EXTERIOR  KNOTS  IN  THE  B-SPLINE  BASIS 
FOR  A  SPLINE  SPACE 


J.  KOZAK 


1.  The  Result. 

In  a  practical  computation  one  is  rarely  able  to  make  statements  about  the  Inverse  of 
a  given  matrix,  particularly  if  the  linear  system  to  be  solved  depends  on  several  free 
parameters.  This  note  is  intended  to  demonstrate  that  in  an  important  linear  problem,  that 
of  spline  interpolation,  the  properties  of  a  chosen  basis  allow  us  to  draw  an  interesting 
conclusion. 

To  start,  let  I  [a,bj  c  R  be  a  given  real  interval,  partitioned  by  the  sequence 


a  «:  t  <  t  <  ...  <  t  <  t  s»  b  , 
k  k+ 1  n  n+ 1 


(1.1) 


with  t^  <  ti+k  ,  all  1  and  some  integer  1  <  k  <  n.  For  the  purpose  of  using  B-splines, 
the  sequence  is  extended  by 


*i  <  <  -  <  Vi  4  fck  '  1 


<  t  <  « •  •  <  t 
n+1  n+2  n+k 

n+k 


(1.2) 


To  simplify  the  distinction  between  both  parts  of  5.  !“  l  Put 

int(£)  s-  (t^)"*^,  ext(t_)  :■  <  t^)  ^  - 1— n+2  '  Th*  colll'ctlon  of  polynomial  splines  of 

order  k  on  I  with  knot  sequence  _t  is  defined  by 

S  ( X )  :■  {f|fM  .  is  a  polynomial  of  degree  <  k  , 

K  I J  * i 1 


jump  f 
i 


(r) 


0,r  <  k  -  card{j|t^  •  t^},  all  i}  . 


Any  f  e  S  (I)  admits  a  unique  B-spline  representation 

K  t 1 

f  -  j,  ai  Ni,k,t 

that  has  proved  very  successful  in  practical  computations.  Here,  (N 


i,k,t 


(1.3) 

)  is  the  B-spline 


partition  of  the  unity,  i.e. 


Ni,k,t(x>  !ti+k  "  ti,fti'ti+1'*“'  ti+k1(*  ”  x>  + 


k-1 
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Fulbright  Grant  No.  79-045-A. 


Consider  now  the  following  interpolation  problem:  for  given  interpolation  conditions 


X  :=  (X^)"=^  and  prescribed  numbers  r  :=  (r^)"=1)  £  Rn  find  f  c  fc(I )  such  that 

X^f  -  r  ,  all  i  .  (1.4) 

There  \  is  an  arbitrary  sequence  of  linear  functionals  with  support  in  I  except  for  the 

assumption  that  the  interpolation  problem  is  correct,  i.e. 

f  e  t<I)  and  X^f  «  0  ,  all  1  ,  implies  f  =  0  . 

The  particular  representation  (1.3)  leads  to  the  solution 

-  (ai I  (1-5) 

with  A  :=  ( X .  N .  ,  ^ ) n  .  ,  .  We  note  that  the  solution  f  does  not  depend  on  ext(t)  , 
t  l  3,k,t  1,3=1  — 

but  A^  and  consequently  a  do.  This  fact  can  influence  the  numerical  accuracy  of  the 
computed  spline. 

Cox  [8]  considered  various  ways  of  choosing  the  additional  knots.  As  he  concluded 
from  numerical  evidence,  in  the  case  of  general-purpose  algorithms  the  coincident  choice  at 
er.d  points  is  preferable.  In  fact,  he  observed  in  [7)  that  the  spectral  condition  number 
of  the  matrix  obtained  in  particular  spline  least-squares  problems  is  considerably  smaller 
for 


ext( t) 


(V  v#' 


t  . ,  t  . 
,  n+ 1  n+ 1 


k-l 


k-1 


(1.6) 


compared  with  an  equidistant  and  an  average  choice.  We  note  also  that  in  the  book  by 
de  Boor  [5]  the  choice  (1.6)  is  the  rule. 

Theorem  1.1.  Let  the  knot  sequences  satisfy 

int(T^)  =  int(t), 

T.  <  t  ,  i  =  1 ,2, •  •  •  ,k-l , 

i  i 

x  >  t  ,  i  =  n+2 ,n+3 , • • • ,n+k  , 

i  l 

and  let  n  >  2(k-1).  If  the  interpolation  problem  (1.4)  is  correct,  then 

n a~  1  ii  >  ii a"1  n  .  d.7) 

T  00  t  00 
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The  inequality  (1.7)  displays  the  numerical  advantage  of  using  the  ext^t)  as  defined  in 
(1.6).  Also,  following  [2],  consider  the  interpolation  map 


Then 


P:  Rn  ♦  S  (I)  : r  +  Pr 


( X  Pr  -  r^  ,  all  i  ) . 


1PI 


sup  (sup|Pr(x)|)  =  sup  ( sup  I  £  a  N.  (x)|)  <  I  a”1  I 
lrl„<1  xel  Irl^O  xel  i-1  1  1,K'-  -  “ 


(1.8) 

and  (1.6)  gives  the  best  bound  of  the  form  (1.8).  It  is  also  clear  from  [2]  that  A  must 

X 

become  singular,  as  or  t  ♦  *.  We  shall  explicitly  observe 


i  .  k-1  t.,,-T  k-1  T  .  -t 

„.-1.  ^  1  r  ~  k+1  r  „  n+r  n  •, 

DA  D  >  -  OcX^  n  -  ,  II  -  J  . 


(1.9) 


I  A. 


t  « 


r=2  “  t  r-2  t_.  -t 

k+1  n+r  n 


Our  final  observation  concerns  complete  even  order  spline  interpolation.  Let  k  »  2m  and 

1  (6  '  ,  •  • « ^,6  .  •  •  * ,  5  ,6^m  ",♦••,6  )  (1.10) 

tk  fck  fck  Vi  tn  Vi  Vi 

with  tf<r,f  f*r\t).  It  is  proved  in  (4)  that  the  Gauss  elimination  without  pivoting  can 
be  applied  safely  if  the  matrix  is  totally  positive.  Unfortunately  for  X  given  in  (1.10) 


the  matrix  is  not  totally  positive.  But  (1.6)  gives  us  A^  of  the  form 


I 


I 


2.  Its  Proof 


Consider  two  knot  sequences  t_  ,  t  that  satisfy  (1.1),  (1.2).  Assume  that  int  (•) 
is  strictly  increasing.  The  general  case  will  follow  from  continuity  properties  of  divided 


differences.  Of  course  there  exists  a  matrix  Q  such  that 

t  2 

A  «  A  Q  , 
r  t  *t  t 


(2.1) 


but  it  is  not  so  obvious  that  we  can  find  Q  ^  explicitly  when  intt^)  “  int(t).  Recall 


Marsden’s  identity  [9] 


(y-x)  -  »  ^lfk#£(y)  Ni(k,t{x)  '  x'y  £I 


Here 


w  =-  (,-Vi>  • 


(2.2) 


(2.3) 


As  observed  in  [1],  (2.2)  implies 


n  k-1 


(y-x 


,k+_1  ’  j,  e-i  <y  "  t4+1>+  *Vk,t(x) 


,  y  e  int(t)  ,  x  e  I 


(2.4) 


But  (2.4)  holds  for  an  y  <  t^  =  a  too,  since  then  both  sides  vanish  identically.  Hence 

lyoyr*’*'yk1  ("  *  x)t_1  “  ?  tyo'y1',,*'yk1  ”  (*  -  VA,k,t(x)  on  1 

i»1  4-1  - 

(2.5) 


under  the  assumption 

k 

Y  !=  *yi*i=0  C  u  1 ^  ‘ 

Similarly,  with  y  c  int(t)  u  ]t  ,•(  , 

—  n+ 1 

(y0'yl',”’yk,(*  ■  x,+_1  =  fyo-yi'”*-yk>(x  -  1 


k-1 


,L  tyo'y,'  — 'V.J  (Vi  '+  -i,k,t'"'  -  -•  (2.6) 


(t„  ,  -  •) .  H,  ,  _ (x)  on  I. 


i=1  "  "  4-1 

One  finds  that  (2.5),  (2.6)  are  the  identity  (3,  (5.10)],  adjusted  to  the  finite  interval 


Lemma  2. 1.  Let  n  >  2 ( k— 1 ) ,  int(T)  =  int(t).  Then 


(Ni,k,t,i-1  ’  CNi ,k  ,t 1 i= 1  Qt  t 


(2.7) 


-5- 


wlth 


1 


Ct  T  !“  T^i.j-I 


and 


T 


[T..,  ”  T.)  ^  V" 


(T  -T)'v"  (t„  )/<?'  (T.  ),  1  t  1  V  1  <  k-1, 

j+k  j  i,k,t  i+k  j-l,k+2,2  1+k 


(T 


5+k 


-  V  £  'i.k.t'V^j-l.k+Z^V  •  n-k+2  ‘  j  <  1  ^ 


i-j 


otherwise. 


Proof.  Assume  for  a  moment  that  int(^)  is  strictly  increasing.  Choose  1  <  j  <  k-1. 

Since  n+1  >  2(k-1)  +  1  »  2k-1  we  can  use  (2.5)  for  any  such  j  to  obtain 
n  k-1 

N  =  l  (T.  -  T.)[x  ,T. T  )  II  (•  -  t  )  N  on  I  , 

3,k,T  ^  3+k  3  3  3+1  j+k  i_1  t+i  +  i,k,t_ 


and  further 


(Ti+k  -  Tj>lTj'Vi',**'Vx1  J!, (*  ■  w+ 


**  X  (T* '  fc^>+ 

‘V*  "  V  J  j+k 

j  ".  <Tr-V 
ft*j 

l+r 


>i«  7,  “r  -  t.*.’ 

0  +  ,Tj+k  -  V  J+k-Fk -  *  qij,t  T  '  aU  1 

n  (Tr  -  V 
*-j 

i*r 


The  proof  for  n-k  +  2  <  j  <  n  follows  in  the  same  way  from  (2.6),  and  (2.7)  obviously  holds 
for  the  remaining  range  of  j  .  The  proof  is  completed  for  a  strictly  increasing 
int(t),  but  int  ( t )  =  int  (t)  simplifies  q  ,  say  for  the  range 

—  —  1  j  T 

1  <  i  <  j  «  k-1  ,  to 


i  j  ,t  i 


k-1 

n 

£=i+1 


( T  j+k  -  V'WWl'-'Vk1'.  n..'*  -  t£,/.n.  '•  -  V’  '  (2-R’ 


k-1 

/  n 

i=i 


and  the  general  case  follows. 


fl 
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2.  Its  Proof 


Consider  two  knot  sequences  t_  ,  x_  that  satisfy  (1.1),  (1.2).  Assume  that  int  (•) 

is  strictly  increasing.  The  general  case  will  follow  from  continuity  properties  of  divided 

differences.  Of  course  there  exists  a  matrix  ^  such  that 

A  =  A  Q  ,  (2.1 ) 

I  t*t  t 

but  it  is  not  so  obvious  that  we  can  find  ^  explicitly  when  int(j)  =  int(t).  Recall 
Marsden's  identity  [9 J 

.k-1 


(y-x)  '  -  ^  k/i(y)  Ni(k(t(x)  '  x'y  £I 


Here 


Vk.t  (*  -  ‘w1  ' 


(2.2) 


(2.3) 


As  observed  in  [1],  (2.2)  implies 


n  k-1 


(y '  x,+-1  ’  J,  <y "  Vi’-  NiA.t(x) 


,  y  £  int(t)  ,  x  e  i 


(2.4) 


But  (2.4)  holds  for  an  y  <  t  =  a  too,  since  then  both  sides  vanish  identically.  Hence 

*  n  k-1 

tyOyr’**'yk,<*  “  X)+”  “  I  [y0'yl'*“'yk'  1 ",  ’  ti+i,+Ni,k,t(X)  on  1 

1=1  £=1  — 


(2.5) 


under  the  assumption 


y  !=  ^yi^  i=o  c  int<i>  u 


Similarly,  with  y  c  int(t)  u  It 

—  n+ 1 


[yO'yl'“''yk1  (*  “  X,+  1  =  ty0'yl',”’yk)  (X  _  *’+  1 


(t.  ,  -•)  N,  ,  ( x )  on  I. 


n  k-1 

”  J,  ly0'yi'“''yk’^J  -Hi  '♦  "i.k.t—  -  -  (2.6) 

One  finds  that  (2.5),  (2.6)  are  the  identity  (3,  (5.10)),  adjusted  to  the  finite  interval 


I. 


Lemma  2.1.  Let  n  >  2 { k— 1 ) ,  int(T)  =  int(t).  Then 


(N,  .  ) n  ,  -  (N,  J?  Q 

i,k,T  1=1  i,k,t  i=1  t  T_ 


(2.7) 


with 


-5 


In  order  to  investigate  Q  further  we  compute,  for  1  <  i  <  j  <  k-1  , 

(Tj+k  '  V  ^i,k,t<t*+k)/'p,j-1,k+2,T  (1 


qij,t  T 


=  (1 


<TUk> 


^  j+k  +  Ti  '  Ti  +  V^i,k,t(tUk,/^j-1,k+2,T(TJt+k> 


'^i,k,t(ti+k)/  *j-1,k+1,T(Tt+k)  +  J1#i,k,t(Vk)/^,k+1,r 


4. 


<W‘ 


and  consequently 
n 


j 


k-1 


=  !  qij<ii=  A  qij,t  x  SLli  v’i,k,t(t4tk)/'(,,k-1,k+1,T(T4+k)  w 


Since  additionally  intt^)  =  int(t)  ,  then 


k-1 

2k-1 

V 

n 

r=1 

L 

Jl-i+k 

2k-1 

n 

r»k 

r*i 

i  r+i 


'l  V 


k-1 

n 

r=  1 


[Tk'Tk+1',”'T2k-11  11  < 


”1-0-1 


k-1 

i+k-1  11  lTi  ‘  tr+i> 

-  t  )  -  f  -JZJ - 

r+i'  2k- 1 


4-k 


n  (T  -  x  ) 
r*k 
r*i 


The  following  lemma  summarizes  this  observation. 

Lemma  2.2.  Let  n  >  2  (k-1),  and  int(x^)  **  int(t)  •  Then 

n 


£  qi),tr=  £  qi3dt  =  1  '  411  1  • 


From  Lemma  2.1  we  conclude 


A  *  K 
T  t  *t  T 


and  to  complete  now  the  proof  of  the  theorem  it  is  enough  to  show  that,  for  the 


-7- 


particular  _t  and  _T  as  specified  there. 


Qt  T  *  0  :=  T  *  0  '  aU  i<j) 


since  then  by  Lemma  2.2  and  invertible  A 


HA-1 1  =  IQl1  a"1  I!  >  — 

T  “  *t  It00  HQ 


t  T 


—  =  II  A~  1  tl 


Also  (1.9J  is  confirmed  just  by  changing  A^  and  ^  in  the  previous  argument.  It  was 

pointed  out  to  me  by  de  Boor  [6]  that  Q  ^  is  known  to  be  totally  positive  if  _t  is  a 

refinement  of  t  on  I  ,  l.e.  i  .t(T)  =  int  (t),  T,  <  t. ,  T  _  >  t  ,  .  The  following  lemma 
—  -  —  k  1  n+2  n+k 

indicates  that  this  fact  holds  in  the  more  general  situation  of  the  Theorem  1.1,  and  the 
columns  of  deserve  to  be  called  discrete  B-splines  as  in  [3]. 

Lemma  2.3.  Under  assumptions  of  the  Theorem  1 . 1 

q(.  „  T  >  0  .  T  „  >  0  ,  all  i.j  .  (2.9) 


Proof .  Choose  again  1  <  i  <  j  <  k-1 .  The  assumption  on  T. ,t  reads 

Ti  <  *i  ‘  *k  <  Ti+k  '  (2-10) 

From  (2.8)  we  observe  that  (2.9)  holds  for  j  =  i  .  Let  j  >  i+1.  The  elements 


q^  ^  are  linear  in  any  t^  ,  i+1  <  i  <  k-1 .  Thus  from  (2.10) 


C  sgn(q 


sgn(q ,  . 


ij,t  T  |t4  - 


ij'l  1  ~  { 


sgn(q 


ij,t  T  |tt 


or  after  repeating  the  argument  k-i-1  times 


SgniqH.t  T>  ’  Sgn,<Jij(t  TJ 


'i+1 


1i+1 

or 


k-1 


m+i 


lk-1 

or 


s9n((Tj+k  ■  ytTi+k'Vk+1'•••'Tj+k1(t•  -  s'  /tnj  (*  -  ,2*m 


for  some  m  ,  j-i-1  <  m  <  k-i-1,  and  some  j<v.  <v.  .<•••<«.  <  k-1  .  For 

J  J  j  j+1  i+m 

m  »  j-i-1  the  argument  of  sgn(*)  in  (2.11)  vanishes,  and  it  is  enough  to  consider 
m  >  j-i  .  Then 

m+i 

tTi+k>  Ti+k+i'“*'Tjy (<’  ■  Vm/  Jj;j  (*  -  \” 

(2-tk>m 

“  tTi+k  Wi**** 'VkVV'V  f'*"’Tv  tK  --z—- 
J  j  3+1  m+i 

(2  -  t  I" 

■  fT  .T  ,•••,!  1  - r 


i+k  i+k+ 1 '  '  j+k  z  ^^m+1+i-j 


,  ,  r  c  m  i,  ■*,  r-m-i+  j-1 ,  *  .  .ra-r 

lTi+k'Wl'*,*'tj+K12  j0  (  r  )<z  ■  w)  (w  *  V 

Hlf-  l-t4 

<j-i)J  rf;0  V  r  <m+i-j-r)l  1  '  '  V 

,  ,j-l(  «•  ,»+i-j  *  .j-1.*  *  ,-m-1 

(“)  (  ^  JU  -  tk)  (w  -  tk)J  (2  -  W) 


for  some  w,z  with  i  <  T  <  w  <  T  <  T,  ,  t  <  2  <  T  .  Thus  the  first 
j  v  k-1  i+k  j+k 

inequality  in  (2.9!  is  confirmed.  The  other  follows  similarly.  ■ 
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